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Abstract 

, Motivated by newly discovered properties of instantons on non-compact spaces, we realised 

' that certain analytic invariants of vector bundles detect fine geometric properties. Wc 

, present numerical algorithms, implemented in Macaulay 2, to compute these invariants. 

' Precisely, we obtain the direct image and first derived functor of the contraction map 

7r: Z — > X, where Z is the total space of a negative bundle over P 1 and tt contracts 
U s ' the zero section. We obtain two numerical invariants of a rank-2 vector bundle E on Z , 

^ ! the width h°(X; (ir*E) vv /ir*E) and the height h°(X; i? 1 ^^), whose sum is the local 

"^h ' holomorphic Euler characteristic x loc (E)- 
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1 Introduction 

In this paper we present effective algorithms, implemented in Macaulay 2, for the computation 
of two numerical invariants of locally free sheaves of rank 2 and with c\ = on a family of open 
complex surfaces which contain a distinguished line £ of self-intersection I 2 = — k, k > 0. 
The interest in these sheaves arises from mathematical physics, since the Kobayashi-Hitchin 
correspondence identifies a certain subset of these sheaves with instantons on Z^, and in this 
picture our two numerical invariants add up to the local charge of the instanton near the line 
I. However, the invariants are strictly finer than the charge, and they apply to a larger class 
of sheaves than just those which correspond to instantons, and they provide a way to stratify 
the moduli of s[2-bundles on Z^ into "nice" components. 

The mathematical theory behind these sheaves and their relation to physics has been 
studied in [GKM08], and the study of their moduli is the subject of [BGK] and [BGK09]. 
The explicit computation of the numerical invariants has been an essential ingredient of 
several of the results in those papers, for the proof of which one used "direct computa- 
tion". It is the aim of this paper to describe general algorithms for these direct computa- 
tions. A reference implementation in Macaulay 2 can be found on the author's website at 
http : //www. maths . ed.ac.uk/~s0571100/Instanton/. 

Finally, we will see that the algorithms are actually easily adaptable to a larger class of 
computations of sheaf cohomology on more general spaces. One such adapted algorithm will 
be used on our upcoming paper [GK], where we study sheaves on local Calabi-Yau threefolds, 
such as Tot(C P i (-1) C P i (-1)) and Tot(C P i O f i (-2)) . Sheaves on Calabi-Yau threefolds 
are of interest to numerous mathematicians and physicists, in particular the aforementioned 
spaces appear in the context of brane theory in papers by Dijkgraaf-Vafa and G. Moore. 



Outline. In §2 we give the basic definitions of the types of sheaves in whose computations 
we are interested in this paper, along with some background results. We define our two basic 
invariants, the width and the height. The algorithms for the explicit computation of the width 
and the height are described in §3 and §4, respectively. Finally, we describe how to compute 
similar invariants of endomorphism bundles in §5 and §6 and how to adapt the algorithms to 
other situations. 



Acknowledgments. The algorithms in this paper are generalisations of the work of I. 
Swanson and E. Gasparim in [GS05]. The author is very grateful to D. Grayson and M. 
Stillman for their great work on creating Macaulay 2 and maintaining a friendly and active 
user community, and to E. Gasparim for the privilege of enjoyable collaboration and sage 
advise. 
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2 Definitions and background 

Let Z k be the total space of the line bundle Opi(— k) over P 1 and k > 0. Denote by I the 
zero section, so that £ 2 = —k. Let £ be a holomorphic rank-2 vector bundle over Z k with 
c\{E) = 0. It is known from [Gas97] that E is an algebraic extension of algebraic line bundles, 

— 0(-j) -^E^ O(j) — , (2.1) 

where 0(j) is the pull-back to Z k of O v i(j) under the projection Z k — > P 1 . 

We fix once and for all coordinate charts U = {z, u} and V = {w, v} on Z k , so that w = z^ 1 
and v = z k u. (A more symmetric picture is given by weighted homogeneous coordinates 
[xq : x\ : X2\ of degrees (1, —k, 1), but this is less useful for explicit computations.) 

In these coordinates, an extension of the form (2.1) is given by the transition function 
from the J7-chart to the T^-chart that takes the form of the matrix 

o p( :~f) ■ w 

The integer j > is called the splitting type of E, and the function p, which is holomorphic 
on U f~l V, represents an element of Ext^ z (0(—j),0(j)). We write \p] for the equivalence 
class of all functions that determine isomorphic extensions, and we say that the pair (j,p) 
determines the bundle E on Zj.. Hence (j, q) defines isomorphic bundles for all q £ [p] . 
It turns out that in our case the function p is always a polynomial: 

Proposition 2.1 ([BGK]). If (j,p) determines a bundle of type (2.1) on Z^, then we can 
choose p to be a polynomial in u, z and z~ x . Moreover, the polynomial can be chosen to have 
the form 

L^J i-i 

P(Z,U) = ^ PrsU T Z S , (2.3) 

r=l s=kr— j+l 

and ifu\p(z, u), the equivalence is given by simply setting all terms that do not appear in (2.3) 
to zero. □ 

Remark 2.2. If u does not divide p, then the extension (2.2) defines a bundle E that is 
also an extension of line bundles of a lower splitting type. In that case the problem (j,p) is 
ill-posed. This does not happen if u divides p. 

We want to compute explicitly from this data the so-called local holomorphic Euler char- 
acteristic 

X loc (E) := h°{X k ; (n*E) vv /ir*E) + h°(X k ; R^E) , 

where is obtained from Z k by contracting the zero-section via tt: Z k — * X k . Since Tr\z k \£ 
is an isomorphism onto X k \{0}, the sheaves (7r*E) vv /ir^E and R^ir^E are supported over 
the single point G X k , and so their spaces of global sections are simply their values at 0. In 
symbols, we have 

H°(X k ; (7r*£) vv /tt*E) = ((7r*£O v 77r*£) =: Q , and 
H Q {R\*E) = {R^E) . 
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To compute these stalks on Xk we make heavy use of the Theorem on Formal Functions and 
instead compute sections of E on Z^. 

Theorem on Formal Functions (Grauert, Grothendieck) . Let ir: Z — > X be a proper map 
of complex spaces and T a coherent sheaf on Z. For x 6 X let £ := tt^ 1 (x). Then 

n 

where denotes the n th infinitesimal neighbourhood of £ in Z. 

3 Computation of h°(X k ; (tt*£) vv /tt*E) 

Let E be a vector bundle on Z^ determined by the data (j,p) as described in the previous 
section. We want to compute the dimension of the vector space Qq := ((7r*£ , ) vv /-7r^£ , ) , 
which is the stalk at of the skyscraper sheaf Q defined by the exact sequence 

— > tt*E (tt,£) vv — > Q — > . 

By definition, the Ox k -module structure on ir*E is determined by the lifting map 

tF: Xk ^n*0 Zk , (3.1) 

which is an isomorphism away from and whose stalk at is, in [/-coordinates, just 7r(u?j) = 
z l u, where 

£>x k ,o =: S = C[w ,wi, . . . ,w k ]/(wiWj - w i+1 Wj-i) , 

where the ideal contains all the indices % = 0, . . . , k — 2 and j = i + 2, . . . , k. 

Thus we are lead to compute the space of sections of E, first as C[z, u] -module and then 
as an S'-module. Since we are using the Theorem on Formal Functions for the computation, 
we will actually be computing the C[[z,tt]]- and 5 A -module structures. However, on any 
Noetherian locally ringed space (X,0) the complet ion O is a flat O-module, and = 
T x <S>e> O for any O-sheaf T and all x & X. 

The computation of tt*E proceeds in three steps: First we apply the Theorem on Formal 
Functions to express (tt*E)q in terms of the cohomology of E, i.e. 

(tt*E)* = \imH°(£^; E^) . 

n 

In local coordinates, elements of H°{£^; E^ n ^ are sections of E of the form a = (a(z, u),b(z, it)) , 
where a and b are a priori power series in 

n 

O ein) (U) = ®u\£{z} . 

i=0 

However, we require that the local section patch correctly onto the other chart, so that 
Tcr = (z>a + pb,z~ib) is a holomorphic section of E^(V), i.e. holomorphic in (z 1 ,z k u). 
This shows that a and b are in fact polynomials, i.e. they contain only finitely many non-zero 
powers of z. 
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Remark 3.1. We have essentially demonstrated the GAGA correspondence for the projective 
schemes £^: all holomorphic sheaves are algebraic; and even if we start with a holomorphic 
section of a sheaf, we are forced to conclude that it is algebraic. Consequently, it is irrelevant 
in our computations of cohomology whether we consider as an algebraic scheme (with the 
Zariski topology) or a complex analytic space (with the Euclidean topology). 

For the second step, we have to show that we can compute the module structure of (it*E))q 
from a finite amount of data (essentially by only going up to a finite infinitesimal neighboor- 
hood, but see §3.2.1). To be slightly more precise, we will not compute H°(^£^; E^), but 
instead we will identify finitely many elements in H°(^£^; E^ n ^ that generate (tt*E))q as 
a C[[z, u]]-module. (The fact that we can do this depends crucially on the structure of the 
space Zfc and the fact that the conormal bundle of £ C is ample.) 

Finally, once we have computed 

M:=H°(£W. £ W) 

(for some sufficiently large N) as a C[[z, u]]-module, the third and final step is to find the 
S^-module structure on M induced by the lifting map (3.1). Here we exploit the fact that u 
is not a zero-divisor in C[[z,u]] and that every element in C[[z, u]] can be expressed in terms 
of Wi = z l u after multiplication by a sufficiently high power of u. 



3.1 An example 

An example is often more illuminating than a detailed theoretic description of an algorithm, 
so let us start with a typical one: 

Example 3.2. Consider Zi, the total space of Opi(— 2), whose blow-down X2 is a surface 
with an ordinary double point, which for convenience we give coordinates x = wq = u, 
y = w\ = zu and w = W2 = z 2 u, where xw = y 2 . Let the E be the bundle on Z2 determined 
by the extension p(z, u) = u and of splitting type j = 3. To compute the module M we are 
thus looking for sections (a, b) of E such that 

( z3 ;r) < 3 - 2 > 

is holomorphic in z" 1 and z 2 u. Since a and b are holomorphic in (2, it), we can write 
a(z, u) = a rs u r z s and b(z,u) = b rs u r z s . 

r,s>0 r,s>0 

The basic idea is to work "one infinitesimal neighbourhood at a time" , i.e. to deal with each 
power of u separately, starting from 0, until one has "enough" information. Thus, starting at 
r = 0, we see from the second entry of (3.2) that z~^b(z, u) (mod u) has to be holomorphic 
in 2 _1 , so that 

b(z, u) = 6 o + bmz + bo 2 z 2 + 603 z 3 (mod u) . 
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Next, there can be no holomorphic terms in a with r = 0. Thus one continues at r = 1: 

a(z, u) = {fliQU + a\\uz + a\ 2 uz 2 + • • • ) (mod u 2 ) . 

The term z 3 a\ou is not holomorphic in z 2 u) , but it is matched by the subsequent ubo^z 3 . 
Thus we have a relation: a±o + &03 = 0. There are no further terms in a(z, u) for r = 1 that 
can be matched by ub(z,u), so ai s = for s > 1. 

We could now carry on to the next formal neighbourhood, find more terms for b 

biou + bnuz + b\ 2 uz 2 + bi3uz 3 + b^uz 4 + bi^uz 5 , 

and then calculate relations on a2 S - However, we shall see immediately that this adds no new 
information to the S'-module of sections. 

We must now find generators of the sections of E, but considered as a module over S, 
where S is the ring 

S = C[x, y, w]/ (xw — y 2 ) . 
We certainly have the following generators: 

O = , Pi = (°) , and fa = ( ° ) . 

But also, we must take into account the relation aio + &03 = 0. Thus the last generator is 

1 = (?) ' 

Of course the module is not free over S, since we have the following relations: 

y(3 = xf3i y(3i = x(3 2 

wf3 = yPi w[3i = y(3 2 

The computation is actually complete now: Even if one were to consider higher generators, e.g. 
the section /?4 = (0, u) coming from the 6io-term, it would just be a multiple of an existing 
generator (here fin = xfio). Also, the a2o-term appears to provide a new, free generator 
a = (u 2 ,0), since z 3 a,2ou 2 is actually holomorphic on V; however, we have a = yfa — x l- 

At this stage of our guiding example we record the S'-module of sections that we just 
computed: 

M := S[p , /3i,&,7] / (yflo - xfi^wfio - yf3i,y(3i - x[3 2 ,w[3i - y(3 2 ) 

We now proceed to compute M v , M vv and the quotient. The simplest, linearly indepen- 
dent elements of M v that we can write down are 

/3 V = (A) !-» x,(3i !->• y,(3 2 (->• w,7 i-> 0} and 7 V = {fa ^ 0, 7 i-> l} . 

A moment's reflection shows that all other possible maps are combinations or S-multiples of 
these two generators, and clearly there are no relations. Thus M v is the free S'-module 

M V = S[/3 V , 7 V ] • 
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The bi-dual is now simply 

M VV = S 
Evaluation on M yields: 



^v = {/5 v^ l}i7 w = {7 v^ l} 



ev(/? ) 
ev(/3 2 ) 



/3 V i-> x 
7 V h-> 

7 V i-f 



x/? vv 



;/3 v -> 



2//3 



vv 



,vv 



So 



coker(ev) = (/? vv ) { 



which has dimension 1. So = 1. □ 
3.2 Description of the algorithm 

The example of §3.1 suggests a general algorithm: We must consider two polynomials a and 
b, use the condition that z^a(z, u) +p(z, u)b(z, u) and z~^b(z, u) be holomorphic in (z~ 1 , z k u) 
to obtain relations on the coefficients a rs and b rs , thence create the S- module M, and finally 
compute the dimension of the quotient M vv /M. 

The crucial consideration is that we only need consider finitely many terms in a(z, u) 
and b(z,u), and this will suffice to describe the module structure of M. In other words, we 
guarantee that we can choose a priori polynomials 



B 



(Xfg % Hi 



and b(z, u) 



r=0 s=0 



B 2 

EE 6 -- 

r=0 s=0 



;Z U 



in which we treat the coefficients a rs and b rs as indeterminates, which toghether with the 
finitely many relations among them generate the module M. The bounds A\, A2, B\ and B2 
will only depend on k, j and p, and they will be determined at the start of the algorithm. 
This is described in §3.2.1. 

Following the computation of the relations among the coefficients, we require a small, 
technical routine to convert the C[z, u]-module into an 5-module. These technical algorithms 
are described at the end of §3.3. Finally, for the computation of the quotient M vv /M we use 
the same computational method that was described in [GS05, Lemma 2.1 (in)]. 



3.2.1 Computation of M 

Intuitively, it is clear that to compute M one has to write down "enough" terms of a and b, 
calculate / := z 3 a + pb and set to zero all terms in / that are not holomorphic in z^ 1 and 
z k u. This gives a set of relations among the coefficients a rs and 6 rs , which in turn determines 
a set of sections that generate M. (In the example of §3.1, the relation 010 + 603 = implied 
the generator 7 = (—it, z 3 ) T .) In this section we give precise instructions on how to find the 
relations among coefficients and how to build from them a generating set of sections. 
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First let us fix some notation: To each coefficient a rs and b rs , let us associate, respectively, 
"elementary" sections 

/ z s u r \ 

a(a rs ) := a® s := I J and cr(b rs ) := a rs : = 

Then the generator associated to a relation R = a rs + Rifiu = 0, where Ru is non-zero for 
at least one (i, I), is a(R) := — tf s + Ylii Ril a u- We denote by 1Z the set of all such relations, so 
we may consider 1Z to be the "solution set" of the holomorphy condition T(" J /T(E; V) = 0. 
With this notation, M is generated as a C[z, u]-module by the set Gn := {cr(R) : R € TZ}, 
and as an 5-module by G'^ := {7r*o"(i?) : R £ 

There are two problems one faces when restricting oneself to a (finite) polynomial, which 
we turn into 




Objectives for the algorithm. 

1. One must find all generators of M, i.e. one must ensure that G-ji generates M. For exam- 
ple, on Z2 with p = and j = 4, the O2o-term contributes a free generator (u 2 , 0), which 
one could miss by only considering the r = and r = 1 infinitesimal neighbourhoods 
for a. 

2. One must find all relations between b r > s i- and a rs -terms. Some 6 rs -terms may appear to 
be free when one does not consider enough a rs -terms. For example, on Z2 with j = 5 
and p = u 2 , the term 6052; 5 may erroneously seem to constitute the free generator (0, z 5 ) 
if one does not include the second infinitesimal neighbourhood and finds 020 + &05 = 0, 
so that the actual generator is (— u 2 , z 5 ). 

There exists a precise bound on the number of infinitesimal neighbourhoods which one needs 
to consider. By including terms from a higher neighbourhoods into the polynomials a or b, 
one may see new relations involving terms from lower neighbourhoods appear, but at the same 
time this will add new generating terms for which one might in turn be tempted to find new 
relations in even higher neighbourhoods. However, we have a priori bounds on the terms in 
a and b that ensure that we compute the correct module structure on M. 

3. It is acceptable for 1Z to contain too many relations involving terms in a. This happens 
when there are not enough terms in b to match. In [GS05] this was called a "fake 
relation" . However, if R £ 1Z is such a fake relation, and if by considering higher terms 
we would find the corresponding "real" relation to be R', then we can ensure that cr(R') 
is already contained in the module generated by G^- 

This will inevitably be the case when p contains several terms of different degree in u: 
In that case one cannot possibly find all correct relations among a finite set of terms. 
The key is to allow high terms of a to be set to zero "erroneously" , rather than to miss a 
relation between a term b r i s > and a term a rs . (The latter would cause us to add a wrong 
generator, while the former only removes a potential generator - but we are careful to 
miss only multiples of earlier generators.) 
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We illustrate this important point by means of Example 3.2: Suppose we only considered 
b up to the neighbourhood r = 0, and a up to r = 2. Then we had to conclude the 
relation R : 022 = 0, which is a "fake relation", whose corresponding "real" relation is 
R' : d22 + bi5 = 0. However, cr(R') = uz 2 a{pf) = ^2(7(7), so we do not need the generator 
a(R'). 

The range of coefficients which one needs to consider depends on the extension p: 
Definition 3.3. Let p € C[z, z~ x ,v\. We define: 

• min n := the minimal degree of u occurring in p, 

• max u := the maximal degree of u occurring in p, 

• min^ := the minimal degree of z occurring in p, and 

• max z := the maximal degree of z occurring in p. 

• If p = 0, then all the above values would be —00; however, for this case we define 
min u := 0, which will later save us from having to consider this case separately. 

For a given bundle E on Z\~ determined by (j,p), there are immediate bounds on the 
number of degrees of z that need to be considered for each fixed r: 

Proposition 3.4. For any degree r of u and independent of p, only the terms 

b r0 u r H h b rM+j u r z kr+J 

occur in b. 

Proof. We require that z~^b(z, u) be holomorphic in z~ l and z k u. By multiplying b r s zS u r 
by z~i we see that the only terms that are holomorphic in z~ l and z k u are those claimed. □ 

Proposition 3.5. For any r < min U; the terms in a of degree r in u, if any, are 

a r ou r + • • • + a T ^ r -jU r z kr ~^ . 

Proof. Since r < min u , no term a rs u r z s can be combined with any term in pb, so the problem 
reduces to making holomorphic in z k u, which results precisely in those terms stated. 

□ 

Proposition 3.6. For r > min^, the only terms a rs u r z s that can possibly be non-zero satisfy 

< s < max{ k(r — min u ) + j + max{max z , 0}, kr — j\ . 

Proof. Consider all terms in z^a rs u r z s that are not holomorphic in z~ l and z k u: They must 
vanish unless they can be matched by a term in pb. The only terms in pb that have degree r 
in u are of the form b r ' s >u r z s , where r — max. u < r' < r — mm u . Since the terms in b are as in 
Proposition 3.4, s has to run at least up to kr' max + j = k(r — min u ) + j, but the multiplication 
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pb may have shifted the term matching a rs by up to max(0,max z ) places up, which explains 
the first term in the statement. 

Secondly, terms up to s = kr — j are automatically holomorphic in the expression z 3 a, so 
if kr — j is greater than the previous expression, all terms up to kr — j must be considered, 
and all the coefficients are free. □ 

Finally, we must turn the Objectives 1-3 into ranges for r that we choose to consider. 
Proposition 3.7. By considering only a truncated generic section 

\ ) - E E ««(%")+ E E - 7 + E E (£) . (") 

u>/iere 

:{ |j'/£f|,max u } + mm u , 

'■-{kr — j, k(r — min u ) + j + max{max 2 , 0}} , and 
:{[j'//c],max u } , 

one finds enough generators to compute the S -module M. 

Remark 3.8. This statement contains two facts: First, we claim that our choice of polynomi- 
als a and b gives enough coefficients from which we form the generators Gn of M. Secondly, 
we claim that the set TZ of relations is correct in the following sense: If we set A = B = oo 
and if 7£oo denotes the associated set of relations, then one of two things happens for each 
R £ TZoo'- Either R is already a relation in TZ, or o~(R) is an S'-multiple of cr(R') for some 
R' G TZ. (This case was illustrated in Objective 3.) 

Proof. Let us denote the three big sums in (3.3) by Ei, S2 and S3 respectively from left to 
right. Since the section (a, b) is holomorphic on U, we must have s > 0, and the upper bounds 
for s in each of the three sums is given respectively by Propositions 3.5, 3.6 and 3.4. To justify 
the choice of the remaining bounds, consider the condition 

a\ _ ( V a + p(z, u) b^ 



a 


:= max 


(3 


:= max 


1 


:= max 



bj-y z ->b 

The term Si is seen to contribute free generators of M, since no term in TEi can be matched 
by any term from TE3. 

The important choice is that of the bound a. Once this has been chosen, we will only 
consider b up to u-degree a — min U) such that p(z,u)b is matched with a (this justifies 
Objective 3 above). This will justify the choice of 7 = a — min u . Moreover this ensures that 
there cannot be any generators coming from o that are erroneously considered as free. It 
remains to prove that our choice of the bound a leads to correct computation of the module 
M. 

By construction, all the generators we get from a are correct, while the generators coming 
from b are either correct or fake. We have to show two things: (i) all fake relations are 
multiples of genuine relations, and (ii) any relation of M is a multiple of a relation that we 
have already found. But both (i) and (ii) follow directly from the choice of a. □ 
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3.2.2 Computation of M vv and l(Q) 

This last section is merely included for completeness. It is no computational obstacle to 
compute the dual and bi-dual of M: 

M v := Hom s (M, S) and M vv := Hom 5 (Af v , S) . 

The evaluation map ev : M M vv is the natural map given by 

ev(a) : 4> h-> 0(a) for all a £ M, (j) £ M v . 

Lastly, note that dimension is invariant under completion, i.e. dirnQg = dimQo, so we have 
/(Q) dim(coker(ev)) . 

3.3 Implementation of the algorithm 

Our reference implementation of the algorithm is written in Macaulay 2 [M2], a computer 
algebra package for commutative algebra. The technical aspects of this implementation are 
specific to that language, and the Macaulay ^-code is available from the author's website at 
http://www.maths.ed.ac.uk/~s0571100/Instanton/. Here we present the generic part of 
the algorithm in pseudo-code. 

Input and output. The algorithm takes as input the data (k,p,j), where k > and j > 
are integers and p is a polynomial in (z ,u). The main function iWidth(k,p, j) computes 
the width l(Q) for the bundle E on determined by (j,p). 

Auxiliary functions. The main function iWidth(k,p, j) calls several auxiliary functions: 
The function makeSectionsAndRing(k,p, j) creates the polynomials a(z,u) and b(z,u) ac- 
cording to Propositions 3.4, 3.5 and 3.7. The function getRelations(k,fTv) computes the 
relations among the coefficients of a and 6, where fTv = z J a(z,u) + p(z , u)b(z , u) . (Note 
that fTv contains all the necessary information.) The function makeModule constructs the 
(S-module M from the data aPoly and bPoly, which arise respectively from a(z, u) and b(z, u) 
by applying all the relations. (For example, if 020 + &05 = is a relation, then we substitute 
Q>20 &05 i n a -) Finally, qLength(M) computes l(Q) from the module M\ for its implementa- 
tion we refer to [GS05]. 

The main function. Name: iWidth. Input: (k,p,j). Output: the instanton width l(Q). 
Pseudo code. 

{aPoly, bPoly, allVars} := makeSectionsAndRing(k, p, j) 
fTv := z"j * aPoly + p * bPoly 
relRes := getRelations (k, fTv) 

apply substitutions from relRes to aPoly and bPoly 
M := makeModule (k, aPoly, bPoly, allVars) 

return qLength(M) 
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Auxiliary function. Name: makeSectionsAndRing. Input: (k,p,j). Output: the 
polynomials a(z,u) and b(z,u), and allVars, a collection of all coefficients occurring in 
a(z, u) and b(z, u). 
Pseudo code. 

minU : = minimal u-degree of p 

maxU := maximal u-degree of p 

minZ : = minimal z-degree of p 

maxZ : = maximal z-degree of p 

aMax := max(ceiling(j/k) , maxU) + minU 
bMax := aMax - minU 

if p = then ( minU = 0; bMax = 0; aMax = ceiling (j /k) ) 

generate coefficients: 
a rs such that r = 0, . . . , minU — 1 and s = 0, . . . , kr — j; 
a rs such that r = minU, . . . , aMax and 
s = 0, . . . , maxjfcr — j, k(r — minU) + j + maxjmaxZ, 0}}; 
b rs such that r = 0, . . . , bMax and s = 0, . . . , kr + j. 

aPoly := sum(a_(r , s) z"s u~r) 
bPoly := sum(b_(r , s) z"s u~r) 

allVars := collection of all coefficients a_(r,s) and b_(r,s) 
return { aPoly, bPoly, allVars } 



Main algorithm. Name: getRelations. Input: (k,fTv) , where 
fTv = z J a(z,u) + p(z,u)b(z,u). Output: A collection of relations like 
{a2o + b 5, «3i + hA + boe}- 

Synopsis. Each relation is the coefficient of a monomial z s u r in fTv for which s > kr. 
Pseudo code. 
rels := { } 

expSet := set of exponents (r, s) appearing in fTv 

for each (r,s) in expSet do 
if (s <= k * r) then continue 
term := the z s u r -term in fTv 

rel := the coefficient of term, scaled to be monic 
rels := rels + { rel } 
end for 

return rels 
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Auxiliary function. Name: makeModule. Input: (k, aPoly, bPoly, allVars). Here 

aPoly and bPoly are the results of substituting all the relations into the original a(z, u) and 

b(z,u). Output: the 5-module M (e.g. its presentation matrix over S). 

Synopsis. Iterating over each coefficient in allVars, we set this coefficient to 1 and all 

others to to get a section (a, b) of E. By multiplying with a high power of u (called uexp), 

we can express (u N a, u N b) as a section of -k^E, and those sections generate M. 

Pseudo code. 



S := makeRing(k) 

Smodule := image (g) : S 1 -> S 2 

N := the maximum of s — kr over all monomial terms z s u r in aPoly and bPoly 

uexp := ceiling(N/k) 

aPoly := aPoly * u~uexp 

bPoly := bPoly * u~uexp 



for each coefficient c in allVars do 

a := aPoly with c = 1 and all other coefficients = 
b := aPoly with c = 1 and all other coefficients = 
Smodule := Smodule + image (^"S'S): S 1 -> S 2 

° \ piStar(b,S) / 

end do 

return Smodule 

This function calls two further auxiliary functions, makeRing(k) and piStar. The first 
one, makeRing(k), returns the quotient ring S ■= C[wq, . . . , Wk]/(wiWj — Wi+iWj-i) for 
% = 0, . . . , k — 2 and j = i + 2, . . . , k. The second function, pi Star, converts monomials 
u r z s into monomials JJi W T i n &i where ^ rii = r and = s. This is possible because 

we multiplied every term by the sufficiently high power u uexp . 

Auxiliary function. Name: piStar. Input: (p,S), where p is some polynomial in u and 
z in which each term is of sufficiently high degree in u, and S is the target ring. Output: 
The polynomial p expressed in ^-coordinates, where w% = z % u. 
Pseudo code. 

res := // this will store the result 

k : = the number k if S is C[u;o, • • • , Wk]/ (uiiWj — Wi + iWj-\) 

//we have variables w_0, . . . , w_k 

for each term t in p do 
degU : = li-degree of t 
degZ : = z-degree of t 
fctr := 1 

if degZ > k * degU then 

error: this term is not convertible! 
end if 
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diff := k 

while (diff != 0) do 

fctr := fctr * w_dif f " (degZ/dif f ) // degZ/diff is integer division 

degU := degU - (degZ/diff) 

degZ := degZ modulo diff 

diff := diff - 1 
end do 

fctr := fctr * w_0~degU 

res := res + fctr * (coefficient of t) 

return res 



At last, we need to compute the length of the module Q, which equals the dimension of 
M vv I 'M as a C-vector space. The computation is performed by the function qLength using 
a presentation matrix for M; the actual algorithm is precisely the one described in [GS05, 
Lemma 2.1 (hi)]. 

4 Computation of h°(Xk m , R l n*E) 

Let E be a bundle on of type (2.1) determined by (j,p). The sheaf R x ti*E is supported at 
the origin, since tt is an isomorphism everywhere else. Therefore H Q \X^; R 1 t^^E) = (i? 1 7r*£') . 
The Theorem of Formal Function gives 

(R^E)^ = limH 1 (£^;E^) . 

n 

However, this limit stabilises at a hnite n, and so we may simply compute the finite-dimensional 
vector space H 1 (Z&; E); then its dimension is the height of E. 

In this section we present an algorithm that produces a basis for H^^Z^; E^ . For this we 
use the Cech description 

1( , T(E;UnV) 
U ^ k ^> ~ T(E;U)®T{E-V) ' 

so we are looking for sections of E on the overlap U (IV modulo sections on either U or V. 
We recall that U and V are affine, and we may consider our sheaves either as analytic sheaves 
over complex spaces or as sheaves over algebraic schemes; both points of view give the same 
results. 

Proposition 4.1 ([BGK, Lemma 2.9]). Let E be determined by (j,p). Then every 1-cocycle 
in H 1 (Zk'i E) can be represented locally over U as 



E E (°S 

r=0 s=kr—j+l 



(4.1) 
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The idea is the following: The vector space H 1 (Z^; Ej is certainly spanned by all the 
monomial cocycles c rs := (a rs z s u r ,0) T from Equation (4.1), so we need to identify which 
linear combinations of the c rs vanish in cohomology. But c rs vanishes in cohomology precisely 
if there is a function b holomorphic on U such that 

T ^a rs z s u r ^ = ^a rs z s+ iu r + pb\ ^ ^ 

is holomorphic on V . (Here T is the transition matrix for E from Equation (2.2).) Since p is a 
polynomial, only finitely many terms in b need to be considered, and we obtain an algorithm. 

First note that if p = 0, then there can be no relations among the c rs , and H (Z^; E) = 
({crs})c- 

Proposition 4.2. If p ^ 0, let mm u be the smallest degree of u appearing in p. To obtain 
H 1 (Zf.; E) , it suffices to check Equation (4.2) for polynomials of the form 

L^atJ " min " kr 

b(z,u) = ^2 b rs z s u r . 

r=0 s=—j 

Proof. This is immediate from the form of p in Proposition 2.1 and Equation (4.2). □ 



4.1 Description of the algorithm 

The algorithm itself consists of two parts: The first part computes all the linear relations 
between the generators c rs ; it returns a list of all basis elements for H^iZ^E} and a set of 
relations, which may contain lots of redundant information. 

The second part of the algorithm takes these sets of generators and relations and reduces 
them to a minimal set of generators and relations. From this new data, we compute the 
dimension of -ff 1 (Z^; E) as the minimal number of generators minus the minimal number of 
relations. 



First part: rinding relations 

1. Let b be as in Proposition 4.2, treating all the coefficients as indeterminates. 

and s = kr — j + 



.7-2 
k 



2. For each monomial z s u r for (r, s) in {(r, s) : r = 0,.. 
1, . . . , — 1|, do the following: 

(a) Let S be the set of all terms in pb with degree (r, s) in (u, z). 

(b) If 5 = 0, then c rs ■= (a rs z s u r ,0) T is an independent generator of H l (Zk] E) . 

(c) Otherwise, if S is non-empty, let B be the set of all coefficients bu appearing in S, 
and let b' = l£B buz l u' Note that z s u r is proportional to at least one term of 
pb' by construction. 

(d) Remove from pb' all terms that are proportional to z s u r and all terms which are 
holomorphic on U; call the result q. 
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(e) Finally, let Q be the set of terms in q that is not holomorphic on V. If Q = 0, 
then the cycle c rs vanishes in cohomology, otherwise we keep c rs as a non-trivial 
generator and obtain the relation z s u r + J2teQ t = 0- 

Implementation. Name: iHeight. Input: (k,p,j) corresponding to the bundle E deter- 
mined by (j,p) on Zj,. Output: a pair (G,R), where G is a set of monomials t such that the 
cycles (t, 0) T span a vector space V, and R is a set of linear relations on V (involving the 
coefficients b rs ) such that H 1 (Zj e ;E) = V / R. 
Pseudo code. 

minU : = minimal u-degree occurring in p 

if p = then minU = 

bMax := floor ((j-2)/k) - minU 

make all indices b_(r , s) for r = 0, . . . , bMax and s = — j, . . . , k * r 

bPoly := ^ rs b_(r,s) * u~r * z~s 

pb := bPoly * p 

aList := list of terms u~r*z~s for r = 0, . . . ,floor((j — 2)/k) and 
s = k * r — j + 1, . . . , —1 

aNonTrivials := {} // These two variables 

aRelations := {} // store the final result 

for each aCycle in aList do 

// To begin, find all terms in p*b that cancel aCycle 

pbpruned := all terms from pb with the same (z, u)-degree as aCycle 

if (pbpruned = {}) then // nothing can cancel aCycle 

aNonTrivials := aNonTrivials + {aCycle} 

continue 
else 

leftb := all terms in bPoly that contain coefficients in pbpruned 



leftpb 
aCycleR 
leftpb 
leftpb 



= leftb * p 

= (all terms in leftpb that are proportional to aCycle) 
= leftpb - aCycleR 

= leftpb - (all terms holomorphic in (z,u)) 



leftovers := those terms of z~j * leftpb that are not holomorphic in (z , z u) 

if leftovers = then 

aNonTrivials := aNonTrivials + {aCycle} 

aRelations := aRelations + {aCycleR + leftovers} 
end if 



end if 
end for 
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return (aNonTrivials , aRelations) 

Second part: reducing to minimal generators and relations The first part of the 
algorithm produces two sets of data: a set G of generating monomials of for form z s u r (i.e. the 
cocycle (z s u r ,0) T is non-trivial in H l (Zk; E^j , and a set R of relations which are polynomials 
in z,z ,u with coefficients bij. Let C be the set of all coefficients bij that can appear; C is 
determined by Proposition 4.2. To find minimal generators and relations, proceed as follows: 

• Build a new set R" of relations without indeterminates as follows: For each relation 
r € R, for each (3 E C, set ft = 1 and all other coefficients in C \ {/?} to zero; add the 
relation r|^ =ljC <\{^} = o to R" . 

• Build a new set of generators G' and a new set of relations R' by starting with G' = G 
and R' = R" as follows: Let N be the set of monomial relations in R', i.e. relations 
of the form z s u r = 0. For each r £ N, remove r from G' and substitute r = into 
every relation in R' . Let iV be the new set of monomial relations in R' and repeat until 
N = 0. 

• The final set G' is a minimal set of generators, and the final set R' is a minimal set of 
relations. 

Implementation in pseudo code. Name: f ixHeightRelations. Input: (G,R), the sets 
of generators and relations which the iHeight algorithm produced. Output: a new pair 
(G' , R' ), where G' is a minimal set of generators for the vector space ^(Z^E), and R' is 
a new set of linear relations, usually empty. Thus |G'| — |R'| is the actual value of the height 
of E (this number is also returned in the actual implementation). 
Pseudo code. 

if (G = {} or R = {» then 
return (G, R) 

end if 

rels := {> // this stores the result R 
allvars := the set of coefficients b_(r,s) 

for each term t in G do 
for each v in allvars do 

11 := t with v=l and all other variables set to zero 
if 11 != then rels = rels + {11} 
end for 
end for 

prunednontrivs := G 
prunedrels := rels 

nullguys := the set of one-term relations (e.g. z 3 u 2 = 0) in prunedrels 
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while (nullguys != {}) do 

for each term t in nullguys do 
replace (R + t) by (R) in prunedrels 
replace (R + t) by (R) in prunednontrivs 
end for 

nullguys := the set of one-term relations in prunedrels 



end while 



return (prunednontrivs, prunedrels) 



5 Note on computing H 1 (Z^; End E^j 



In the next two sections we compute invariants of the endomorphism bundle £nd(E) = E<S>E A . 
This bundle plays a fundamental role in the deformation theory of the sheaf E: H 1 (Z).;£nd E) 
is precisely the tangent space at E of the moduli of holomorphic bundles diffeomorphic to E. 
This follows for example from [DK90, Proposition 6.4.3], as the Kuranishi map vanishes on 
Z k (since H^Z^T) = for alH > 2 and every coherent sheaf J-) . 

If determines a bundle E on Z\. as before, with transition function T given by 

Equation (2.2), then the endomorphism bundle End E = E (g> E v is a rank-4 bundle whose 
transition function is given, after a convenient change of coordinates 
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The computation of H 1 {Z\.\£nd E) is essentially the same as for H 1 {Z^\ E) , in the sense that 
there exists a simple canonical representative for every 1-cocycle, and the first part of the 
algorithm from §4 translates almost literally, while the second part remains unchanged. An 
implementation of this algorithm, named hlend(k,p, j), is contained in our Macaulay 2 code 
at http : //www. maths . ed.ac.uk/~s0571100/Instanton/. 



6 Cancelling infinities: the computation of H°(Zk;£ndE) 

On the other hand, the space PP^iZ^EndE^ is infinite-dimensional, so we cannot compute 
it directly. However, if we are given two different bundles (j,pi) and (j,P2), we can compute 
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a "relative dimension" of H -spaces as follows: For each n > 0, the space 
is finite-dimensional, so we can compute the difference A n (pi,p2) '■= h° (^ n ); E(pi)( n )} — 
hP [l^ 1 ; E(p2)^) • Since p±, pi are polynomials, A n is constant for n S> 0. Finally, we can 
define a function 

A(E)=h(j,p) := lim A n (0,p) , 

which is non-negative since the split bundle given by p = has the largest amount of sections 
on each infinitesimal neighbourhood & l \ (See [BGK] for a discussion of this non-trivial fact.) 

In the remainder of this section we describe an agorithm to compute h° (l^ ; End E 1 ^) from 
the input data (j,p) and n. This amounts to finding the most general section a € Y\£nd E;U) 
such that Sa € T(£nd E;V) . Let a = (a,b,c,d) T , where we write a typical component as 
a(z,u) = J2 r s>o a rsZ s u r . We have 



S 



b 
c 
\dj 



I a + z° pb + z ipc + p 2 d\ 
z 2 ^b + zipd 



-■2.1 



c+ z 3 pd 
d 



(6.1) 



For computing H°, we see that d and we know that z J p is holomorphic on V, thus so is 
J pd. This implies immediately that 



kr 



d(z, u) 



^2 d rs z s u r 

r>0 s=0 



kr+2j 

and c(z, u) = c rs z s u r 

r>0 s=0 



Next, in the second entry of (6.1), z 3 pd contains powers z s u r for s < kr + 2j — 1, so that we 
may take 6(z, u) = X^r>o Ss^y 1 ( we could have taken into account the fact that u divides 
p for an even sharper bound, but choose not to). Finally, similar considerations show that 
we may take a(z,u) = J2 r >o Ss=o 2j_1 o-rsZ s u r . Here we assume that j > and that p is in 
canonical form (2.3); if j = 0, there exists only one bundle anyway. 

We now describe an algorithm that computes for a given n > the vector space 
The algorithm generates a set of linear relations on a larger vector space spanned by mono- 
mials, and in our reference implementation we use a built-in function from Macaulay 2 to 
compute the dimension of the resulting space directly. 



Implementation. Name: hOend. Input: (k,p, j ,n), where (k,p,j) determine as before 
a bundle E on and n is the infinitesimal neig hbourhood. Output: h°(i^;E^). 
Pseudo code. 

aPoly := Er=o ES)^ a -( r . s ) z ~ s u ~ r 

bPoly := E™=oEs=o b - (r ' s) z ~ s u " r 

cPoly := Er=oEs=o 2j c -( r > s ) z ~ s u ~ r 

dPoly := E™=oEs=o d - (r ' s) z ~ s u " r 
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el := aPoly + z~j * p * bPoly + z~(-j) * p * cPoly + p"2 * dPoly 

e2 := z~(2*j) * bPoly + z ~j * p * dPoly 

e3 := z~(-2*j) * cPoly + z"(-j) * p * dPoly 

e4 := dPoly 

relations := {} 

for each polynomial pol in (el, e2, e3, e4} do 
s := z-degree(t) 
r := u-degree(t) 
for each term t in pol do 
if s <= k * r then 

continue 
end if 

badterms : = terms in pol of the same degree as t 
badcoefs := badterms / z"s * u~r 
relations := relations + {badcoefs = 0} 
end for 

return (the dimension of (a rs , b rs , c rs , d rs ) /relations) 

7 Note on adapting the algorithms to other spaces 

To conclude this paper we outline how to adapt the algorithms from this paper to any space 
X = Tot(0™ =1 Opi(oi)) for which we have a GAGA principle and in which we can contract 
the zero section. We refer to [BGK] and infer that X satisfies GAGA at least when we have 
a,i < for all i, in which case all vector bundles are algebraically filtered, and in particular 
every rank-2 bundle is an algebraic extension of algebraic line bundles. 

The changes in the algorithm are basically as follows: We now need coordinate charts 

U = {(z,v,!, . . . ,u n ) and V = {(w,vi, . . . ,v n )} 

which patch together via w = z -1 and Vi = z~ a *Ui. If E is a bundle on X of type 
(2.1), then its transition function is determined by the splitting type j and a polynomial 
p € C[^ ±:L , tti, . . . , u n ], which can be put into a canonical form. 

Example. On the Calabi-Yau threefold W\ := Tot(C P i(— 1) ® C P i(-l)), the polynomial p 
has the canonical form 

2j-2 2j-2-t j-l 

p(z,u,v) = Y^ ^2 Yl Ptr S z s u\u\, e€{l,2}. 

t=e r=l—e s=r+t— j+1 

The putative local sections of E and End E that one uses in the computation of l(Q) 
and K (X; End E) are now also polynomials in z, u±, . . . , u n , but otherwise the algorithms are 
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essentially identical. The adapted algorithms for the flop space W\ are available together with 
the previous algorithms on the author's website. 

We wish to state a last result, on which the author hit only after adapting the algorithm 
to W\ and always computing zero. 

Proposition 7.1 (cf. [GK]). Let X be as above, with contractible zero section isomorphic to 
W , and let tt: X — > X' be the contraction. If diva X > 2, then for any locally free sheaf £ on 
X, the width of £ vanishes, i.e. 

H°(X';(n*£) vv /(■**£)) =0 . 

A Application: Rank-2 Bundles on and genericity 

In this section we tabulate some results that were computed with the algorithms from this 
paper. We will consider bundles of rank 2 on the surfaces Z^ := Tot(Opi(— k)) of splitting 
type j which are determined by a polynomial p according to Equation (2.2). The moduli of 
such bundles has been studied in [BGK], and here we provide concrete examples of elements 
of the moduli and their numerical invariants. 

This section is divided into two cases, first the instanton case for which j = nk for some 
positive integer n, followed by non-instanton bundles. This terminology arises because it 
was shown in [GKM08] and [BGK] that a holomorphic bundle of the type which we consider 
corresponds to an instanton on Z^ precisely when j = nk (via a version of the Kobayashi- 
Hitchin correspondence). 

A.l The instanton case 

In the case that j = nk for some n € N, the bundle E determined by (j,p) corresponds to 
an instanton on Z^. In this case it was shown that the split bundle is the unique one with 
the highest invariants (h,w) (the height and width), while the generic bundles are those with 
the lowest invariants (h,w). Our computations indicate that the condition of being generic 
or split also corresponds to h l {£nd E) being respectively minimal or maximal. 
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In fact, the numerical evidence leads us to conjecture one and discover another relation 
between the invariants (h,w) and (h 1 , A). 
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Conjecture: w + h = \ = ((/i 1 — A) — j)/2 + j/ k. 
Proposition A.l. A + h 1 = h 1 (End (split)") , or equivalently 

h l (£ndE) - h° {End E\ l(m) ) = h 1 (End(split)) - h° (End (split] £{m) )) 
for all large m. 

Proof. We can express the statement in terms of the Hilbert polynomial 

(f>E,m(n) := x{£nd E {m) (n)) = h l (End E\ i{m) (nj) - h°(£nd E\ i(m) (n)) ; 

then the statement is 4>E,m(fy = 0spiit,m(O). But in fact we have <pE,m(n) = 0spiit,m("-) for 
any n and m, since the Hilbert polynomial cannot distinguish extensions on £( m ' C Zj., as 
one sees from direct computation. We retain the clause "for all large m" since h 1 (End E\^m)) 
stabilises to h 1 (End E\ eventually. □ 

From [BGK] and the Conjecture we get h 1 + A = n(2nk + k — 2) for j = nk. We also get 
that for the generic bundle p = zu, we have h 1 — A = (3/e — 2)n — 2 for k > 1 and h 1 — A = n 
for k = 1. 

Corollary: Assuming the Conjecture, we have A = n 2 k — x, h 1 = kn(n + 1) — 2n + x- 

A. 2 Non-instanton bundles 

Here j ' ^ (mod k). 
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3 


4 


u + Z 3 U 


10 


2 


8 


(0,3) 


3 


4 


z 3 u 2 
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z~ x u 
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zu, Z 2 U 


14 
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10 


(0,4) 



Continued on next page 



23 



Computations of instanton invariants 



A: 
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h 1 -A 


(w,h) 
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We see that the invariants (w, h) are now required to distinguish the generic bundles (those 
with lowest (w,h), wheras the split bundle is no longer the only one with the highest values 
of (w,h) (see [BGK] for details). However, h}\£ndE\ still distinguishes the split bundle. The 
physical interpretation of these "non-instanton" bundles invites further exploration. 

B Usage example 

The algorithms that we described in this paper are implemented in Macaulay 2 version 1.1, 
the code is contained in the file InstantonInvariants2 .m2. Suppose we want to study the 
bundle on of splitting type 7 given by the polynomial p = z~ l u + zu 2 . We set Macaulay 
up as follows: 

$ M2 InstantonInvariants2 .m2 
Macaulay 2, version 1.1 
with packages: . . . 

11 : p = z*u~2+z~-l*u; 

Now we compute the width and height, and h l {Zi\ £nd E^j . When calling iWidth, we use 
the option Verbose=>f alse to suppress additional output. 

12 : iWidth(2,p,7,Verbose=>f alse) 

13 : iHeight(2,p,7) 

14 : f ixHeightRelations oo 

15 : hlEnd(2,p,7) 

16 : f ixHeightRelations oo 

We find that the width is 2, the height 6, and h 1 [Z2; £nd E) = 33. 
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